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Abstract 

The probability density function of velocity fluctuations of glanulence observed 
by Radjai and Roux in their two-dimensional simulation of a slow granular flow 
under homogeneous quasistatic shearing is studied by the multifractal analysis for 
fluid turbulence proposed by the present authors. It is shown that the system of 
granulence and of turbulence have indeed common scaling characteristics. 
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1 Introduction 

In this paper, we apply the multifractal analysis (MFA) [1, 2, 3, 4] of fluid turbulence to 
granular turbulence (granulence [5]) in order to see how far MFA works in the study of 
the data observed by Radjai and Roux [5] in their two-dimensional simulation of a slow 
granular flow subject to homogeneous quasistatic shearing. Radjai and Roux reported that 
there is an evident analogy between the scaling features of turbulence and of granulence 
in spite of the fundamentally different origins of fluctuations in these systems. MFA is 
a unified self-consistent approach for the systems with large deviations, which has been 
constructed based on the Tsallis-type distribution function [6] that provides an extremum 
of the extensive Reny [7] or the non-extensive Tsallis entropy [6, 8] under appropriate 
constraints. 
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2 Multifractal Analysis 



MFA of turbulence rests on the scale invariance of the Navier-Stokes equation for high 
Reynolds number, and on the assumption that the singularities due to the invariance 
distribute themselves multifractally in physical space. 

The velocity fluctuation Su n = \u(»+£ n )—u(») \ of the nth multifractal step satisfies the 
scaling law \u n \ = \Su n /Su \ = 5 r ^ 3 with S n = £ n /£o — ( n — 0, 1, 2, • • •). We call n the 
multifractal depth which can be real number in the analysis of experimental data. We will 
put 5 = 2 in the following in this paper that is consistent with the energy cascade model. 
At each step of the cascade, say at the nth step, eddies break up into two pieces producing 
the energy cascade with the energy-transfer rate e n that represents the rate of transfer 
of energy per unit mass from eddies with diameter £ n to those with £ n+ ±. Then, we see 
that the velocity derivative \u'\ = lim^oo with the nth velocity difference u' n = 5u n /£ n 
for the characteristic length £ n diverges for a < 3. The real quantity a is introduced in 
the scale transformation [9, 10] x — > x' = \x, u — > u' = \ a ^ 3 u, t — > t' — A 1_a / 3 £, p — > 
p' = \ 2a / 3 p that leaves the Navier-Stokes equation du/dt + (u ■ V)u = —Vp + vV 2 u of 
incompressible fluid invariant for a large Reynolds number Re = 8ui n £ m /u. Here, v is the 
kinematic viscosity, p = p/p with the thermodynamical pressure p and the mass density 
p, and 5u in and £ in represent, respectively, the rotating velocity and the diameter of the 
largest eddies in turbulence. The largest size of eddies is, for example, about the order of 
mesh size of the grid, inserted in a laminar flow, which produces turbulence downstream. 

Within MFA, it is assumed that the singularities due to the scale invariance distribute 
themselves, multifractally, in physical space with the Tsallis-type distribution function, 
i.e., the probability P (jl \a)da to find in real space a singularity with the strength a within 
the range a ~ a + da is given by [2, 11, 12] P^ n \a) = (ZW) _1 {1 - [(a - a ) /Aa} 2 } 71 /^ 
with (Aa) 2 = 2X/[(1 — q) In 2]. Here, q is the entropy index introduced in the definitions 
of the Renyi and the Tsallis entropies. This distribution function provides us with the 
multifractal spectrum f(a) = 1 + (1 — q)~ l log 2 [l — (a — a ) 2 /(Aa) 2 ] which, then, produces 
the mass exponent 

r(q) = l-a q + 2Xq 2 (l + y^)" 1 + (1 - ^[l - log 2 (l + yjc 9 )] (1) 

with Cq = 1 + 2g 2 (l — g)Xln2. The multifractal spectrum and the mass exponent are 
related with each other through the Legendre transformation [10]: f(a) = aq + r(q) with 
a = —dr(q)/dq and q = df(a)/da. 
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The formula of the probability density function (PDF) (u n ) of velocity fluctuations 



is assumed to consists of two parts, i.e., II^l 



u, 
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where the first 



term is related to P^ n \a) by lTg (\u n \)du n oc P^ n \a)da with the transformation of the 
variables \u n \ = 5 n a ^ 3 , and the second term is responsible to the contributions coming 
from the dissipative term in the Navier-Stokes equation violating the invariance under 
the scale transformation given above. Then, we have the velocity structure function 
in the form {(\u n \ m )) = J du n \u n \ m U( n \u n ) = 2 7 &) + (1 - 2^)a m 5 n ^ with 2 7 £> = 
/ du n \u n \ rn ATl( n \u n ), a m = {2/[y / C m/3 (l + ^/C m/3 )]} 1/2 and the scaling exponent Cm = 
1 — r(m/3) given with the mass exponent (1). 

The PDF n( n )(£ n ) both of velocity fluctuations and of velocity derivative to be com- 
pared with observed data is the one defined through f[( n \£ n )d{; n = Yl( n \u n )du n with 
the variable £ n = u n / {{u 2 ^} 1 / 2 scaled by the standard deviation of velocity fluctuations. 
For the velocity fluctuations larger than the order of its standard deviation, £* < 
(equivalently, |a| < a*), the PDF is given by [3, 4] 
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with £ n , = CnC /3 ~ C2/2 and IlW = 3(1 - 2^ l> ) / {2£ n ^j2nX\ \n8 n \). This tail part rep- 
resents the large deviations, and manifests itself the multifractal distribution of the sin- 
gularities due to the scale invariance of the Navier-Stokes equation when its dissipative 
term can be neglected. The entropy index q should be unique once a turbulent system 
with a certain Reynolds number is settled. For smaller velocity fluctuations, |£ n | < 
(equivalently, a* < \a\), we assume the Tsallis-type PDF of the form [3, 4] 



U, 



n (ri) <! i - ^ (i + 3/V)) 



dU r 
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(3) 



where a new entropy index q' is introduced as an adjustable parameter. This center part 
is responsible to smaller fluctuations, compared with its standard deviation, due to the 
dissipative term violating the scale invariance. The entropy index q' can be dependent on 
the distance of two measuring points. 

The two parts of the PDF, (2) and (3), are connected at £* = £„5"* /3 ~ C2/2 with the 
conditions that they have a common value and that their slopes coincide. The value a* 
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is the smaller solution of C2/2 — a/3 + 1 — f(a) = 0. The point £* has the characteristics 
that the dependence of n*" - -*^*) on n is minimum for n 1. With the help of the second 
equality in (3) and (2), we obtain AU.^(x n ), and have the analytical formula to evaluate 
7^. Their explicit analytical formulae and the definition of £ n are found in [3, 4]. 

3 Turbulence 




Figure 1: Analyses of the PDF's of velocity fluctuations (closed circles) and of velocity derivatives 
(open circles) measured in the DNS by Gotoh et al. at R\ = 380 with the help of the present the- 
oretical PDF's Z7"( n )(£ n ) for velocity fluctuations (solid lines) and for velocity derivatives (dashed 
line) are plotted on (a) log and (b) linear scales. The DNS data points are symmetrized by 
taking averages of the left and the right hand sides data. The measuring distances, r/rj = i n /r/, 
for the PDF of velocity fluctuations are, from the second top to bottom: 2.38, 4.76, 9.52, 19.0, 
38.1, 76.2, 152, 305, 609, 1220. For the theoretical PDF's of velocity fluctuations, \i = 0.240 
(q = 0.391), from the second top to bottom: (n, n, q') = (20.7, 14.6, 1.60), (19.2, 13.1, 1.60), 
(16.2, 10.1, 1.58), (13.6, 7.54, 1.50), (11.5, 5.44, 1.45), (9.80, 3.74, 1.40), (9.00, 2.94, 1.35), (7.90, 
1.84, 1.30), (7.00, 0.94, 1.25), (6.10, 0.04, 1.20), and £ = 1.10 ~ 1.43 (a* = 1.07). For the the- 
oretical PDF of velocity derivatives, (n, h, q') = (22.4, 16.3, 1.55), and £* = 1.06 (a* = 1.07). 
For better visibility, each PDF is shifted by —1 unit in (a) and by —0.1 in (b) along the vertical 
axis. 

The dependence of the parameters «o, X and q on the intermittency exponent [i is 
determined, self-consistently, with the help of the three independent equations, i.e., the 
energy conservation: (e n /e) = 1 (equivalently, r(l) = 0), the definition of the intermit- 
tency exponent /i: (e 2 n /e 2 ) = (equivalently, ji — 1 + r(2)), and the scaling relation: 
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1/(1 — q) = \joL- — l/at+ with a± satisfying f(a±) = 0. Here, e is the energy input rate 
to the largest eddies. The average (• • •) is taken with P^ n \a). 

The PDF's extracted by Gotoh et al. from their DNS data [13] at R x = 380 are 
shown, on log and linear scales, in Fig. 1 both for velocity fluctuations and for velocity 
derivatives, and are analyzed by the theoretical formulae (2) and (3) for PDF's. We found 
the value \i = 0.240 by analyzing the measured scaling exponents ( m of velocity structure 
function with the formula given above, which leads to the values q = 0.391, a = 1.14 
and X = 0.285. Through the analyses of the PDF's for velocity fluctuations in Fig. 3, we 
extracted quite a few information of the system [14, 15, 16, 3, 4]. Among them, we only 
quote here the dependence of q' on r/rj: q' = — 0.05 log 2 (r/?y) + 1.71 [4]. 

4 Granulence 

Let us now analyze the velocity fluctuations in glanulence simulated by Radjai and Roux 
[5]. Since they observed that the fluctuations share the scaling characteristics of fluid tur- 
bulence, we try to investigate the system by means of MFA which extracted, successfully, 
the rich information out of turbulence as was seen in the previous section. The power 




Figure 2: Analysis of the experimental PDF of fluctuating velocities, measured in the quasistatic 
flow of granular media by Radjai and Roux, with the help of the present theoretical PDF Z7"( n )(£ n ) 
for velocity fluctuations (solid lines) are plotted on (a) log and (b) linear scales. The experimental 
data points are symmetrized by taking averages of the left and the right hand sides data. The 
integration time r, normalized by a shear rate, for the experimental PDF are, from top to 
bottom, 10~ 3 , 10 _1 . For the theoretical PDF, u = 1.347 (q = 0.930), from top to bottom: 
(n, q') = (88.0, 1.28), (40.0, 1.22), and £* = 1.14, 1.14 (a* = 0.364). For better visibility, each 
PDF is shifted by —1 unit in (a) and by —0.1 in (b) along the vertical axis. 
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spectrum of the fluctuating velocity field on one- dimensional cross sections exhibits a clear 
power-law shape with the slope —f3 with (5 ~ 1.24 [5], which is quite similar to the power- 
law behavior with the slope —5/3 in the inertial range of the Kolmogorov spectrum [17]. 
However, the granular model is an assembly of frictional disks, the power-law observed in 
granulence does not mean the energy conservation in contrast with the case of the energy 
cascade model for fluid turbulence. 

For the conditions to determine the parameters cto, X and q, we adopt, instead of 
the energy conservation, the slope of the power spectrum, i.e., /3 = 1+C 2 = 2 — t(2/3) 
in addition to the definition of the intermittency exponent and the scaling relation. The 
latter two are the same as those for turbulence. As there is no experimental data, for the 
present, to determine the intermittency exponent fi for granulence, we cannot have the 
values of the three parameters through the three conditions. Therefore, we determine the 
value of the intermittency exponent by adjusting the observed PDF with the theoretical 
formulae (2) and (3), since the accuracy of the formulae in the analysis of PDF's for 
turbulence is quite high as was shown in the previous section. The best fit of the observed 
PDF of fluctuating velocities by the formulae (2) and (3) is shown in Fig. 4. We found 
the value [i = 1.347 giving q = 0.930, a = 0.377 and X = 0.050. By making use of the 
mass exponent with these values, we have (e n /e) = S~ T ^ with r(l) = 0.648 representing 
a breakdown of energy conservation. It is attractive to see that the result is quite close to 
(e n /e) = 3/2 which may be consistent with the coefficient of friction 0.5 for the simulation 
[5]. We further extract the relation between r and t n as r = 1.3 5„ ' 131 by comparing the 
observed flatness and the one with the theoretical PDF's (2) and (3). This relation may 
be a manifestation of the fact that Taylor's frozen turbulence hypothesis does not work 
for granulence. 

5 Prospects 

We showed with the help of MFA that the system of turbulence and of granulence have, 
actually, common scaling feature in their velocity fluctuations as was pointed out by Radjai 
and Roux [5] . We expect that various observation of granulence will be reported at higher 
statistics, and that one can extract more information out of the data to determine the 
underlying dynamics for granulence in the near future. 
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